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Abstract 

K'^S'^T [5] recently derived a new 6*^-order wave equation i('(il/6; {d'^ + '&Uxdx+^Uxx){'U't + 
"^xxx ~\~ 6?^^.) = 0, found a linear problem and an auto-Bfickclund transformation for it, and con- 
jectured its integrability in the usual sense. We prove this conjecture by constructing an infinite 
commuting hierarchy KdVn^ with a common infinite set of conserved densities. A general con- 
struction is presented applicable to any bi-Hamiltonian system (such as all standard Lax equations, 
continuous and discrete) providing a nonholonomic perturbation of it. This perturbation is conjec- 
tured to preserve integrability. That conjecture is verified in a few representative cases: the classical 
long-wave equations, the Toda lattice (both continuous and discrete), and the Euler top. 



1 Introduction 

The theory of differential and difference Lax equations has been well understood by the 
middle of the 1980s, and no surprises have disturbed the contented tranquility of the 
subject ever since. Until now. 

Recently, the 5 authors of [5] subjected to the Painleve analysis the 6*'^-order nonlinear 
wave equation 

Uxxxxxx ~\~ OjUxUxxxx ~\~ bUxx'^xxx ~\~ CU^Uxx~\~ 

+duu + eUxxxt + fUxUxt + gutUxx = 0, (1.1) 

where a, b, c, d, e, / and g are arbitrary parameters, and they have found 4 cases that pass 
the Painleve test. Three of these were previously known, but the 4*'' one turned out to be 
new (eqn (5) in [5]): 

[dl + 8uxdx + 4.Uxx){ut + Uxxx + Qui) = 0. (1.2) 

This equation, as it stands, doesn't belong to any recognizable theory. The K'^S'^T convert 
it, in the variables v = Ux, w = ut + Uxxx + 6u^, into (eqn (12) in [5]): 

Vt + Vxxx + l2vVx - tfa; = 0, Wxxx + ^VWx + 'iwVx = 0. (1.3) 

If integrable, this is a truly remarkable system: since {w = 0} leaves only the unperturbed 
KdV itself, and the constrain on w is differential, what we have here is a nonholonomic 
deformation of the KdV (= KdV2) equation. The 5 authors of [5] found a linear problem 
and an auto-Backlund transformation for their equation, but reported that they were 
unable to find higher symmetries on available computers, and asked if higher conserved 
densities and a Hamiltonian formalism exist for their equation. 
All of these queries are resolved affirmatively below. 
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We now proceed to the general construction of nonholonomic perturbations of bi- 
Hamiltonian systems. Rescaling v and t in the equation (1.3), we get: 

Uf = 6uUx + Uxxx — Wx, Wxxx + ^UWx + 2UxW = 0. (1-4) 

This can be converted into 

«t = Si(^)-BiH= B^(^-^)-B\w), B\w) = 0, (1.5) 
where 

B^ = d = dx, B'^ = d^ + 2{ud + du) (1.6) 
are the two standard Hamiltonian operators of the KdV hierarchy, n = 2, and 

Hi =u, H2 = H3 = nV3 - ul/2, ... (1.7) 
are the conserved densities. 



And that's it. The ansatz (1.5) provides a nonholonomic deformation of any bi-Hamiltonian 
system. The question, naturally, is whether this ansatz is reasonable or an absurd phan- 
tasy. My answer is two-fold: (A) It is reasonable; (B) It is difficult, if not impossible, to 
prove integrability in general. The arguments are as follows. 

(A) Each system (1.5) has an infinite sequence of H^'s as its conserved densities: 



dt 5u dt 5u 5u Su 



= B\^^^)w ~ --^^B\w) = 0, (1.8) 

where, as usual, a ^ 6 means: (a — 6) G Imd (a "trivial Lagrangian" ) . 

(B) The above calculation is about the only one that can be reliably performed in the 
{u\ w;}-picture, because the constraint B^{w) = is nonholonomic. Thus, if we proceed to 
develop the variational calculus in the {u; ?ij}-variables, we would be blocked, because the 
calculus works only when the factor /d{^l^) is a free module, where is the module of 
differential 1-forms (see [8]). Thus, the question of integrability: whether the flows (1.5) 
still commute between themselves, can not be answered in general with the modern tools. 
It can be answered for the KdV case (and I believe for all the standard differential Lax 
equations) through a subterfuge. To get a hint on how to proceed, we start in the next 
Section with the classical long-wave system 

ut = hx + uux, ht = (uh)x- 

Section 3 is devoted to the KdVn^ hierarchy (1.5,6) itself. Section 4 treats the Toda lattice 
and its continuous limit. The last section considers the classical Euler top. 
Remark 1.9. The term nonholonomic is of a recent vintage, and seems to have been 
invented by Hertz, see [3]. 
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2 The Classical Long- Wave Equations 

The classical long wave system is bi-Hamiltonian (in fact, 3-Hamiltonian) [11], [6], [8]: 

u\ (h + u'/2\ f0d\f5/6u\fv?h + f^\ 

(2d du \f5/5u\}uh\ 
-\ud hd + dh ) \6/ShJ \2 )■ ^^-'^^ 

Thus, its perturbation (1.5) is: 

ut = {h + u'^/2)x-W2,x, ht = {uh)x - wi^x, (2.2a) 
{2wi + uw2)x = 0, uwi^x + {hd + dh){w2) = 0. (2.26) 

The first of the constrains in (2.26) is resolvable, but the second one is not, and we seem to 
be stuck. The help comes from the missing from (2.1) (gravity) parameter g [1], rescaled 
away for mathematical simplicity (which was immaterial in the holonomic framework, but 
is fatal in the nonholonomic case): 

u\ Jgh + u''/2\ fO d\ (d/du\ fv^h + gJf\ „ , 

2gd du \ (5/5u\ (uh\ 

ud hd + dh ) \d/5h)\ 2 )' ^ 

so that now 

ut = {gh + /2 - W2)x, ht = {uh-wi)x, (2.4a) 
{2gwi + uw2)x = 0, uwi^x + (hd + dh){w2) = 0. (2.46) 

The constraint (2.46) is resolvable as a regular series in g. This can be seen as follows. 
The first eqn in (2.46) yields: 2gwi + UW2 = function of t and e only, and we rescale that 
function into 1: 

2gwi + UW2 = 1 ^ (2.5) 
W2 = {l-2gwi)lu^ (2.6) 

w,,x = —{hd + dh)l{l - 2gw,) = ( - -^i)^ + 25( + (^,). (2.7) 



Set now 



u;i=^/zfc. (2.8a) 



k=0 
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Then 



ZQ = -h/u'^, Zk+i = 2{^d + d^){zk) => (2.86) 



Wi 



f; ('\+ yih/u^^^ = 1[1 - (1 - 44)-V^] (2.9a) 



fc=o 
1.. . /i 



W2 = -{l-^g^r'/^ => (2.95) 



But G commutes with all the HnS, because 

Ghh/Guu = 29/h = H^,hh/Hn,uu, VnGZ>2. (2.11) 

Thus, all the flows (1.5) commute also. 

The workable approach to our general problem hence is this: 
(A) Rescale the variables u in (1.5) in such a way that the nonholonomic constrain B'^{w) = 
becomes resolvable, hopefully in the form 

w = SG/Su. (2.11) 

{B) If then 

{G,iJ„}~0, Vn, (2.12) 

then all the flows (1.5) commute between themselves. (See the end of Section 4 for more 
on this.) 

Let's see now how this approach works for the KdVn6 case. 
Remark 2.13. Since the long-wave system is three- Hamiltonian, the nonholonomic con- 
struction applies not only to the pair {B^,B^), but also to the pair {B^, B^) of the corre- 
sponding Hamiltonian structures. It's not clear how these two different perturbations are 
related. 

Remark 2.14. A'^-component systems of hydrodynamical type (= 0-dispersion) are triv- 
ial for N < 3, but their honholonomic perturbations are no longer so. 

3 KdVn6 

We rescale dt and by e. The KdVnQ (1-4-7) becomes (1.5), now with 

= d, b'^ = e^d^ + 2{ud + du), (3.1) 
B'^{w) = e^w^^ + 2{ud + du){w) = 0. (3.2) 

To solve (3.2), set 



oo 

W = 

k=0 
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We get: 

Wo = u~^^'^, wi = (3-3) 
and in fact 



^ = ^' ^ = E^'^'"'(^f)^'- {G,iJn}~0, Vn, (3.4) 



where ps arc certain differential polynomials from the differential algebra <Q[u^^\u^'^\ ...] [e]. 
The proof is, unfortunately, rather long, and I omit it ([10]). I believe that similar rescaling 
works for the general differential Lax (= Gel'fand-Dickey) hierarchy, with the Lax operator 



TV 



jC = u + ^Ui{edy, UN = 1, UN-i = 0, iVGZ>3, (3.5) 
1=1 

but I haven't proved it. (The method was explained, for the case of the Burgers hierarchy, 
in my talk at the AMS meeting at Williams College in the fall of 2001.) 
Remark 3.6. The situation becomes much more complicated when one passes to the 
modified Lax equations. For the KdV6 case, with the Miura map u = evx — v^, one gets: 

vt = e^Vxxx - Qv'^Vx + P, (3.7o) 
{ed - 2v){p) = d{w), {ed - 2v)d{ed + 2v){w) = 0, (3.76) 

so that one has a pair of nonholonomic constrains attached to one scalar field v. 

4 The Toda Lattice 

The Toda lattice is a classical mechanical system with the Hamtiltonian 

2 

H = ^{^ + e9"+i-9"). (4.1) 

n 

In the variables 

an = Pn, bn = exp{qn+i - Qn), (4.2) 
the motion equations become: 



t 



(1 - A-i)(6)\ _ / (1 - A-i)6 \ fa\ _ f 



b{A-l){a) J \b{A-l) J \lj \b{A-l) 

_ f bA-A-H a{l- A-i)6 \ fl\ _ f bA - A-H a{l - A'^)b \ fS/Sa 



(1 - A-i)6 \ fS/Sa\ a2 
) \5/5b) 4 ' 



6(A-l)a 6(A- A-i)6 ; V b{A-l)a b{A - A-^)b ) \5 / 5b ) ^"'^ ' 
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(4.3) 

where A is the shift operator: (A/)(n) = f{n + 1), all the equalities are understood 
as between functions of n G Z (or Z/NZ), and (4.3) shows the first two (out of three) 
Hamiltonian structures, in the {a; fe}-variables, of the Toda lattice (see [7]). 
The nonholonomic deformation ansatz (1.5) produces: 

at = (1 - A-i)(6 - bw2), bt = 6(A - l)(a - wi), (4.4) 
(6 A - A-^ b){wi) + a(l - A-^){bw2) = 0, (4.5a) 
6(A - l){awi + (1 + A-^){bw2)) = 0. (4.56) 

The nonholonomic constrain (4.5), as it stands, is unresolvable. To proceed, we first rescale 
b into eb and then look for solutions regular in e. We get: 

(6 A - A-i b){wi) + a(l - A-^){w) = 0, (4.6a) 
b{A-l){awi + e{l + A-^){w)) =0, w:= bw2. (4.66) 

(4.66) implies that awi + e(l + A~'^){w) = function of t and e only, and we rescale it into 
1: 

au;i + e(l + A~^)(tt;) = 1 ^ (4.7) 
wi = ^[1 - e{l + A-^){w)] 

-(1 - A-^){w) = i(6 A - A-i 6)^[1 - e(l + A-^){w)] = 

-e{cA-A-^){l + A-^){w), c:=b/aa^^\ (4.8) 

where 

g(^) := A^(g), s e Z. (4.9) 
Setting 



w = — 

k=o 



J2^ke', (4.10) 



we find: 

(l-A-i)(zfc+i) = (cA-A-ic)(l + A-i)(^fe), zo = c, keZ>o. (4.11) 



The latter equation, being nonlocal, looks rather impenetrable; it's not even clear if it's 
solvable. So let's pass to the continuous limit to see what the situation is in simpler 
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circumstances. The previous formulae become: 

f hr.\ /O dh\ f a\ /O db\ ( bjbu \ ( c? 



bj^ \ha^) \hd Qj\lJ \bd 0j\6/6bJ\2 
bd + db adb \ f bd + db adb \ ( 5 /5a 



+ b]= (4.12) 



hda 2bdb J \oJ \ bda 2bdb J \6/6b, 

at = (e6(l - W2))x, h = b{a - wi)^, (4.13) 

{bd + db){wi) + awx = 0, w:= bw2, (4.14a) 

bd{awi + 2ew) = 0. (4.146) 

(4.146) resolves into 

awi + 2€w = l^ (4.15) 
wi = {1 — 2ew)/a 

-w^ = i{bd + db)-{l - 2,w) = {^d + d^){l - 2ew) = 
a a a^ a'^ 

= Ca;-2e{cd + dc){w), c:=b/a^. (4.16) 
Setting w = — YlV=o ^k^^ again (4.10), we find: 

Zk+i,x = 2e(c5 + dc){zk), zo = c^ (4.17) 

''=cty'^'' (4.18) 

and the calculation identical to that of §2 shows that (now with e = 1) 

wi = ry~^/^ ry := - 46, (4.19a) 

W2 = ^ . (4.196) 



Since 

^ = ^ = 2ry-3/^ (4.20) 



db da 

there exists a Hamiltonian G = G{a, 6) such that 
^wi\ f6/6a 



\o/ob/ 

and this G commutes with all the conserved densities -fern's of the continuous Toda flow 
(4.12), because 

Gaa/mf{G) = ^= Hm,aa/ibdt,f{Hm), Vm. (4.22) 

Thus, the continuous limit picture is manageable. Back to the discrete play, the equa- 
tion (4.11.) There exist no general methods to handle nonlocal recurrencies such as (4.11) 
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save for the method of bi-Hamiltonian systems (see [9].) So, we first move (4.11) into a 
skewsymmetric form, by applying from the left the operator (1+ A), resulting in: 

(A - A-'){zk+i) = (1 + A)(c A - A-i c)(l + A-'){zk), zo = c. (4.23) 

UnfortTinatcly, the form (4.23), as it stands, is not of the bi-Hamiltonian character, be- 
cause, e.g., 

= c^^^c + CC + cc(^) Im{S), (4.24) 

i.e., zi is not SH/6c for any H. This however, shouldn't be the end of the story, and it 
isn't. The help comes from the observation that 

zjc = |-(c2 + c(-^)c + cc(i))/2, (4.25a) 
oc 

D{z2/c) = [D{z2/c)]\ (4.256) 

so that Z2/C Im{S); here D{ ■ ) is the Prechet derivative of ( • ). This strongly suggests 
that we set 

Zk = cpk, k G Z>o, (4.26) 
multiply eqn (4.23) from the left by c, and rewrite (4.23) as 

[c(A - A-i)c](pfe+i) = [c(l + A)(c A - A-i c)(l + A-')c]{pk), po = 1- (4.27) 

Miraculously, and for no discernible reason: 

(a) The matrix (in fact, scalar) 

= c(l + A)(c A - A-i c)(l + A-i)c (4.28) 
is Hamiltonian; 

(b) The pair of Hamiltonian matrices, (4.28) and 

b^ = c(A - A-^)c (4.29) 

form a Hamiltonian pair. The bi-Hamiltonian theory then guarantees the existence of a 
sequence of Hamiltonians {hm} such that 

Pm = 5hm/5c, meZ>o. (4.30) 

Thus, our constrain (4.8) has been resolved: 



w = bw2 = -J2^''^i^^k/Sc), c = b/aa^^\ (4.31) 
fe=o 

Setting 

00 

h = J2^''hk, (4.32) 



A;=0 
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and noticing that 



S _ 5 _ I 5 

6a a (5c' Sb aa(^) dc^ 



(4.33) 



we see that 



.6hk 



1 5h 5h 



1 \ ^ p. ( 

-"^2 = ^2^ce — = = — , 

k=0 

w^i = - - e-(l + -c— = - - — = — % a - h). 
a a oc a da da 



(4.34a) 
(4.346) 



Thus, 



Wl 
W2 



fS/Sa 
[6/6b 



{log a — h), 



(4.35) 



and the last step now is to show that loga — h commutes with all the conserved densities 

i^m's of the full Toda lattice (4.3). But this is true in general: 

(C) Suppose that the constrain B^{w) = is (1.5) has been resolved as 



6G 

w = some G. 

du 



(4.36) 



Then 



6Hrr. 

5u 



6G 
6u 



6K 



m—l 



5u 



(4.37) 



Thus, our prescription for analyzing the nonholonomic deformation, stated at the end 
of Section 2, works perfectly provided the nonholonomic perturbation it; is a variational 
derivative (4.36). 

Remark 4.38. With hindsight, one readily sees that the bi-Hamiltonian pair (4.29) 
and 6^ (4.28) describes the classical Volterra lattice 



c = c(c«-c(-i)), 



(4.39) 



see [7]. 

Remark 4.40. The nonholonomic perturbation (1.5) of the Volterra lattice (4.39), 



c = c(A - A-i)c(l - w), (c A - A-^ c)(l + A-^){w) = 0, 



(4.41) 



is not simplified by the rescaling c i-^ ec, but it does so upon the rescaling c h- > 1 + ec 
around the stationary solution {c = 1} of the Volterra lattice. 
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5 The Euler Top 

The constrain B'^{w) = is, in general, nonholonomic only for systems which are either 
differential or difference on Z, i.e., for dimensions 1 and "1/2". In dimensions, i.e., in 
Classical Mechanics with a finite number of degrees of freedom, the constrain B^{w) = 
becomes holonomic. Thus, e.g., there would be no problem in resolving that constrain for 
the periodic Toda lattice on I^/NTj. 

Let's see how this works in practice, for the simplest possible case, the so(3) Euler top: 

xi = aiX2X3, ±2 = a2XsXi, xs = 03X1X2, (5.1) 

where x = {xi,X2,xs) G M^. I follow the notation and results of a short and concise 
presentation in [12]. Let's agree that {ijk) stands for an even permutation of (123). 
Equations (5.1) can be the rewritten as 

Xi = (XiXjX}~, (^-2) 

with a = (ai, 0:2, 0:3) G being an arbitrary but fixed vector of parameters. If c = 
(ci,C2,C3) is another vector in M^, then 



1 ^ 

Hc = ^Y.'^xi (5.3) 

i=l 

is an integral of (5.1) iff 

(a,c)=0, (5.4) 

so that the space of integrals is two-dimensional. 

Now, for any vector 7 G R^, consider the Poisson brackets 

=7fcXfc. (5.5) 

Then 

{xi, i/c}^^^ = CjXjjkXk - CkXkjjXj = {cjjk - Ck^j)xjXk, (5.6) 

so that the conditions 

a = c X 7 ^ (q, c) = 0, (5.7) 

guarantee that the motion equations (5.1) are Hamiltonian in the Hamiltonian structure 
(5.5) with the Hamiltonian (5.3). So, let's choose vectors /3 and 7 such that {/3,7, 

from a right orthonormal basis. Set and B^ corresponding to /3 and 7, respectively: 

-P3X3 -P2X2 

-B^ = I P3X3 pixi \,-B^= \ 73x3 71X1 I . (5.8) 




-Pixi 
The constrain B^{w) = becomes: 




w X X^^^ = 0, X^'^) := (710:1,72X3,73X3), 



(5.9) 
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so that 

w = constX^'^\ const = constat), (5.10) 
and the perturbed motion equations (1.5) become: 

Xi = aiXjXk - const{X^'^'^ x X^t))^. (5.11) 

But 

(X(^) X xW)i = (/3 X 'y)iXjXk = ^^aiXjXk (5.12) 

Thus, finally, the perturbed top equations are 

Xi = const' aiXjXk, (5.13) 

so that the overall effect of the perturbation amounts to the time rescaling of the orig- 
inal top. This is reminiscent of the general Chaplygin theorem identifying some special 
nonholonomic systems with the time-rescaled Hamiltonian ones (see [4,3].) 
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